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Applying the dynamical coherent potential approximation to a simple model, we have 
systematically studied the carrier states in ^i-^Mn^B-type diluted magnetic semicon- 
ductors (DMS's). The model calculation was performed for three typical cases of DMS's: 
The cases with strong and moderate exchange interactions in the absence of nonmag- 
netic potentials, and the case with strong attractive nonmagnetic potentials in addition 
to moderate exchange interaction. When the exchange interaction is sufficiently strong, 
magnetic impurity bands split from the host band. Carriers in the magnetic impurity 
band mainly stay at magnetic sites, and coupling between the carrier spin and the lo- 
calized spin is very strong. The hopping of the carriers among the magnetic sites causes 
ferromagnetism through a double-exchange (DE)-like mechanism. We have investigated 
the condition for the DE-like mechanism to operate in DMS's. The result reveals that 
the nonmagnetic attractive potential at the magnetic site assists the formation of the 
magnetic impurity band and makes the DE-like mechanism operative by substantially 
enhancing the effect of the exchange interaction. Using conventional parameters we have 
studied the carrier states in Gai-^Mn^As. The result shows that the ferromagnetism is 
caused through the DE-like mechanism by the carriers in the bandtail originating from 
the impurity states. 

KEYWORDS: diluted magnetic semiconductors, exchange interaction, coherent potential approximation, 
carrier-induced ferromagnetism 



§1. Introduction 

For more than two decades diluted magnetic semiconductors (DMS's) have attracted much at- 
tention because of the combination of magnetic and semiconducting properties. In A^_ x Mm x B^ 1 - 
type (II-VI-based) DMS's, Mn impurities substituting for 2+ cations act as stable 2+ ions and 
therefore there are few carriers, making these DMS's insulators. It is widely accepted that in II- 
VI-based DMS's a carrier (s electron or p hole) moves over many sites while interacting with the 
localized (d) spins on Mn sites through the sp-d exchange interaction. ^ The exchange interaction 
strongly enhances the effect of magnetic field on band splitting, leading to spectacular magnetoop- 
tical effects (e.g., giant Faraday rotation or Zeeman splitting). In recent years the attention has 
also been focused on III-V-based DMS's (Gai-^Mn^As and Ini_. r Mn r As) due to the high potential 
for new device applications. It is highly noteworthy that doping of Mn into GaAs and InAs leads 
to ferromagnetism, and interesting magnetooptical and magnetotransport phenomena. This ferro- 
magnetism is generally called "carrier-induced ferromagnetism" because hole carriers introduced by 
Mn incorporation mediate the ferromagnetic coupling between Mn ions. 2 ) The microscopic mech- 
anism for carrier-induced ferromagnetism is still controversial. The following properties, however, 
seem to be generally accepted as for (Ga,Mn)As: (i) Mn ions substitute randomly for Ga cations in 
the zincblende structure. 2 ) (ii) A Mn ion in GaAs gives rise to an acceptor level at about 0.113 eV 
above the valence band. 3 ) (iii) The Mn ion has highly localized d states with a magnetic moment 
of ~ 5/iB (or S = 5/2). 3 ~ 5 ) (iv) The Mn-induced states near the Fermi energy play a key role 
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in the origin of ferromagnetism. According to photoemission studies, 6-8 ) X-ray absorption spec- 
troscopy, 9 ) and band calculations, 10,11 ) those states are mainly created in As 4p orbits, (v) The 
p-d exchange interaction between the As Ap hole and the localized d spin is antiferromagnetic, 12 ' 13 ) 
and its amplitude is not very different from that in II-VI-based DMS's. 6,14 ) (vi) Antisite defects 
(As ions sitting on Ga lattice site), for example, are common in semiconductor samples grown by 
low-temperature molecular beam epitaxy. 15 ) Many holes may be trapped not at Mn acceptors but 
at such defects, though we may expect one hole donated by a Mn atom. The density of the holes 
and that of the Mn ions are therefore regarded as separate sample-dependent quantities that are 
to be determined experimentally. 

Since the fabrication of III-V DMS's much effort has been made on the elucidation of the origin 
of the ferromagnetism. To date, the Ruderman-Kittel-Kasuya-Yosida (RKKY) mechanism, 16 ) the 
double-exchange mechanism based on d electron hopping, 17 ) the spin-wave approximation for the 
kinetic-exchange model, 18 ) the double-resonance mechanism, 19 ) the Zener-model description, 20 ) a 
theory based on the random-phase approximation and the coherent potential approximation, 21 ) 
and the coherent potential approach to a model with randomly distributed Ising spins 22 ) have been 
proposed. To our knowledge, however, no systematic theoretical study on the carrier states in 
DMS's, in which both the exchange interactions and nonmagnetic attractive potentials exist, have 
been accomplished taking the effect of disorder into account. 

In this study we aim to clarify the nature of carrier states in DMS's applying the dynamical 
coherent potential approximation (dynamical CPA) to a simple model. Furthermore, on the ba- 
sis of the numerical calculation, we discuss the origin and the mechanism of the carrier-induced 
ferromagnetism in III-V DMS's. 23 ) 

The main content and the organization of this paper are as follows. In §2, we briefly summarize 
the present model and the dynamical CPA. In §3 after a general consideration of the DMS's 
and the model parameters (§3.1), we first present the results for the carrier states in the DMS's 
with no nonmagnetic attractive potential: the case of the strong exchange interaction (§3.2) and 
the case of the moderate exchange interaction (§3.3). Then, based on the result we discuss the 
mechanism of carrier-induced ferromagnetism in DMS's (§3.4). The result of this study reveals 
that the double-exchange (DE)-like mechanism 24 ' 25 ) is in effect to realize a high Curie temperature 
when the exchange interaction is strong enough to form a magnetic impurity band split from 
the host band. On the other hand, the RKKY mechanism may be relevant when the exchange 
interaction is so weak as not to produce a magnetic impurity band. The latter corresponds to the 
case of II-VI-based DMS's. As mentioned above, the strength of the p-d exchange interaction in 
III-V DMS's is not very different from that in II- VI DMS's. Then, what causes the ferromagnetism 
in III-V DMS ? The key lies in the attractive Coulomb potential between a p hole and a Mn 
acceptor center. In §3.5, we investigate the role of the attractive potential to elucidate the origin 
of the carrier-induced ferromagnetism in III-V DMS's. Section 4 is devoted to (Ga,Mn)As together 
with a comparison with (In,Mn)As. Concluding remarks are presented in §5. 

§2. Basic consideration 

2. 1 Model Hamiltonian for a carrier in a DMS 

In order to study the effect of the p-d exchange interaction between a carrier (p hole) and lo- 
calized magnetic moments together with magnetic and chemical disorder in DMS's, we previously 
proposed a simple model. 23,26 ) In the model, the local potentials of nonmagnetic A ions in a semi- 
conducting compound AB are substituted randomly, with mole fraction x, by the local potentials 
that include the exchange interactions between carrier spins and the localized spin moments on 
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magnetic M ions. Thus, the potential to which a carrier is subject at a site differs depending on 
whether the site is occupied by an A ion or an M ion. The Hamiltonian H is given by 

H = ^ ^ ^mn0^ m ^nn ~\~ ^ ~] U n , (2-1) 

m,n,fi n 

where u n is either or uff , depending on the ion species occupying the n site: 

U n = E A^2 a n/i a nti , (2-2) 



Em a nM a ™M ~~ I X/ a nfJ. tJ /*" ' Sndnu ■ (2.3) 



The notations here are conventional and the same as in the previous papers. 23 ' 26 ~ 28 ) Here, Ea 
(Em) represents a nonmagnetic local potential at an A (M) site. We employ a local potential 
Em although the real potential for a hole is an attractive screened Coulomb potential exerted by 
a Mn 2+ acceptor center in III-V-based DMS's such as Gai-^Mn^As. The exchange interaction 
between a carrier spin and the localized spin S n of the Mn site n is expressed by — lander ^ u ■ S n a nu , 
where represents the element of the Pauli spin matrices. 

2.2 Dynamical coherent potential approach 

A carrier moving in a DMS described by Eq. (2.1) is subject to disordered potentials which 
arise not only from substitutional disorder but also from thermal fluctuations of d spins through 
the p-d exchange interaction. Furthermore, when magnetization arises, the effective potential for 
the carrier differs according to the orientation of the carrier's spin. In the dynamical coherent 
potential approximation (dynamical CPA), 29 ' 30 ) the disordered potential is taken into account in 
terms of the spin-dependent effective medium where a carrier is subject to a coherent potential Ej 
or Ej according to the orientation of its spin. The coherent potential E-j- (EjJ is an energy (cj) 
dependent complex potential. Then, a carrier moving in this effective medium is described by the 
unperturbed Hamiltonian K: 

K = J> fc + r>^. (2.4) 

Thus, the perturbation term V(= H — K) is written as a sum over each lattice site: 

V = J> n , (2.5) 

n 

where v n is either v£ or vff , depending on the ion species occupying the n site: 

vi = ^(£ A -S>t M a nM , (2.6) 

v n = ^2( E M ~ E^a^a^ - I^4 M o> • S n a nu . (2.7) 
Next, using the reference Green's function Go given by 

(n%n \G (u>) \nv) = (^m^ ^ ^ K a \i^j > ( 2 - 8 ) 
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where (O)o denotes the expectation value of O in the vacuum, we define the t-matrix t A which 
represents the multiple scattering of carriers due to an A ion embedded in the effective medium by 



ti = (2.9) 



and t M due to an M ion by 



# = (2-10) 



We have omitted suffices of the matrices G, t and v since they are self-evident. The matrix t£ 
ifn) represents the scattering by an isolated potential v£ {vff) in the effective medium completely. 
Note that K, and thus Go, includes no localized spin operator. According to the multiple-scattering 

theory, the total scattering operator T of a random medium, which is related to G = (a a) 

\ u-H , Q 

as 

G = G + G TG , (2.11) 

is expressed as a series, 

n n m (-L n ) n m (^ n ) I (-L m ) 

In the CPA, the condition 

(in)av = for any site n (2-13) 

determines the coherent potential =t,|) and we approximate (G) av by Go- Here we express 

the average of a quantity O over the disorder in the system as (0) av . Since the present system 
includes both substitutional disorder and the thermal fluctuations of the localized spin an M site, 
the average of the t matrix is written as 

(tn)av = (l-*)tf + *<#(S)>. (2.14) 

Here (1 — x) and x are the mole fractions of A and M atoms, respectively; (t M (S)) means the 
thermal average of t M over fluctuating localized spin S. The coherent potential is decided such 
that the effective scattering of a carrier at the chosen site embedded in the effective medium is zero 
on average. Since the spin off-diagonal element becomes zero after the average, the condition (the 
dynamical CPA condition) is given by 

(1 -x)tf T +x(tff(5*)> = 0, (2.15a) 
(1 - s)tft + z(tJf(Sg> = 0. (2.15b) 

For simplicity, the t matrix elements in the site representation < n[i\t\nv > [n is a site index, 
fi, v =| or I) are written as t^ v . It is worth noting that in the calculation of (i^) the spin 
flip processes are properly taken into account, and a single ^-matrix element t£f depends on both 
£f and Ej. Therefore we need to solve Eqs. (2.15a) and (2.15b) simultaneously. Note that the 
diagonal matrix element involves S z as an operator, where S z stands for the z-component of 
the localized spin. The thermal average over the fluctuating localized spin is taken as 

s 
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where h denotes the effective field felt by the localized spins. Since there is one-to-one correspon- 
dence between (S z ) and the parameter A = h/ksT, we can describe the carrier states in terms of (S z ) 
instead of A. Note that the thermal average of off-diagonal t-matrix elements (t*f ) av = (tH) aY = 
because the magnetization is assumed to be along the z-axis. In this work we treat the localized 
spins classically for simplicity. The actual calculations were performed for S = 400. 

Throughout this work we assume a model density of states with a semicircular profile whose 
half-bandwidth is A, 

as the unperturbed density of states. Then the density of states with \i spin D^{u)) is calculated as 
DJcj) = -— Im(ra/i|G(a;)|n//) a v = --Im(n/x|G (w)|n/x) 

TT TT 

= --lm[ A de ^— - (2.18) 

by using S M determined by the condition (2.15). In all of the present results we confirmed numer- 
ically the sum rule: 

/OO /'OO 
D T (w)dw = / D^duj = 1 . (2.19) 
-oo J —oo 

2.3 The Curie temperature 

We calculate the Curie temperature T c to investigate the condition for the occurrence of fer- 
romagnetism. Throughout this paper we assume that the carriers are degenerate. Then we obtain 
the density of the carrier with fx spin n M and the total energy E((S Z )) as 

n M = f F D^)duj (2.20) 

J — oo 

and 

E((S Z }) = r coiD^ + D^du , (2.21) 

J — oo 

respectively, as functions of the Fermi level ep- Here the dependence on (S z ) is contained in D^{uj). 
Note that E((S Z )) is the sum of the kinetic and the exchange energies. For a fixed value of {S z )/S, 
the total carrier density n (= + rij) has one-to-one correspondence with ef and therefore E({S Z }) 
can be expressed as a function of n. The free energy per site of the system is given as 

F({S Z )) = E((S Z ))-TS, (2.22) 

where the entropy due to the localized spins is given by 

s 

S = xk B \og ^2 exp (XS Z ) - xk B \{S z ) . (2.23) 
s z =-s 

The parameter A is determined so as to minimize F((S Z )) through the condition -^F((S Z )) = 
0. If we expand F((S Z )) in terms of ((S z )) 2 , T c is determined as the temperature where the 
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coefficient of ((S z )) 2 vanishes. As is revealed later [see Fig. 20(a)] the energy difference between 
the ferromagnetic and paramagnetic states, 

E((S Z )) — E(0), is approximately proportional to ((S z )) 2 up to full polarization at a fixed 
carrier density. Therefore we fitted the numerical data of [E((S Z }) — E(0)]/A to the expansion 

wm-wd = _j(|>y +6 f<|>v (2 , 24) 



(2.25) 



and estimated the coefficients a and b. Finally we obtain T c as 

k B T c 2a 
A 3x ' 

§3. Carrier states in DMS's 

3.1 General consideration 

Let us start with confirming that the exchange interaction term, — la ■ S, has two energy 
eigenvalues according to the manner of coupling between the carrier spin and the localized spin. 
The parallel coupling state (denoted by p) with 25 + 2-fold degeneracy has the energy eigenvalue 
e p = —IS, while the antiparallel coupling state (denoted by a) with 25-fold degeneracy has the 
energy eigenvalue e a = I(S + 1). The classical spin approximation assumes 5 to be infinite while 
keeping IS constant. Then the parallel and antiparallel eigenstates have the same degeneracy and 
their energy eigenvalues have the same absolute value \IS\. In this work we adopt the classical 
spin approximation since the magnitude of the localized spin on a Mn 2+ (5 = |) is pretty large. 
It is, therefore, sufficient to appoint the value of the exchange energy IS instead of appointing the 
values of I and 5 separately. 

In most DMS's, the p-d exchange interaction between a p hole and a localized spin favors 
antiparallel coupling. The p-d exchange interaction is known to play a crucial role in magnetooptical 
effects in II-VI-DMS's, and is also believed to cause the carrier- induced ferromagnetism in III-V- 
DMS's. Hence, keeping the p-d exchange interaction in mind, we assume IS < hereafter. The 
present model is characterized by only two parameters, IS/ A and Em /A, since we take Ea as the 
origin (= 0) of the energy. 

3.2 Result for IS = —A and Em = 

In the following we show the results for three typical cases of Ai- x M x B-type DMS's, in which 
5 % of the nonmagnetic (A) ions are randomly substituted by magnetic (M) ions. The present 
and next subsections reports typical cases in Em = 0. In Figs. 1 ~ 4 we show the numerical 
results for IS = —A and Em = 0.0, which is the case where the exchange interaction is so strong 
that magnetic impurity bands appear. In Fig. 1, the spin-polarized DOSs, D^(u>) and Z?j_(a>), are 
depicted for various values of (S z )/S. In the dilute limit (x — ► 0) impurity levels appear at the 
energies of ^ = ^ e mtIS ^ _j_ l ^_A_^ _ ±1.25. When x = 0.05, impurity bands form around 
the impurity levels. The total number of states of each impurity band is x, irrespectively of {S z ). 
The low (high) energy impurity band corresponds to the antiparallel (parallel) coupling state. The 
impurity bands are strongly affected by the change in (S z ). On the other hand, the host band 
is negligibly affected. In order to elucidate the origin of the change in the DOS, we calculate the 
species-resolved DOS. In Fig. 2, we depict the A- and M-site components of the DOS, (1 — x)D^(u>) 
and xDjf(u), where D^(uj) [Djf(ui)] represents the local DOS with ^ spin (// =f or j) associated 
with the A (M) ion. The total number of A- and M-site states are 1 — x and x, respectively. Since 
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D^(uj) and Djf(u) are normalized 

D,(u) = (l-x)Di(u;) + xDff(u). (3.1) 

The result shown in Fig. 2 reveals that the impurity state is mainly composed of the M-site states 
and that the change in the impurity band is mainly ascribed to the change in Djf(oj). 

Next, we investigate the manner of coupling between the carrier spin and the localized spin. In 
Fig. 3(a) we show the result for the optical carrier spin polarization P(u>), where P(u) is defined 
by 

P{UJ) ~ dMTd^) ■ (3 - 2) 

The result strongly suggests that 

P(cu) & + ^ ^ av in the lower impurity band, (3.3a) 

P(u>) — — in the higher impurity band. (3.3b) 

We can deduce the result of Eq. (3.3) by assuming that the carrier spin always couples antiparallel 
(parallel) to the localized spin in the lower (higher) impurity band. To confirm the above picture, 
we calculate the spin-coupling strength Q(uj) defined by (see Appendix) 



s _(S(u-H)a.S)/S 



(6(co - H)) 



(3.4) 

M-site 



Thus, Q(uj) corresponds to — (cos 6), where 9 is the angle between the carrier spin and localized 
spin at an M site. The result for Q(u) shown in Fig. 3(b) clearly indicates that 

Q(oj) ~ 1.0 in the lower impurity band, (3.5a) 
Q(oj) ~ —1.0 in the higher impurity band. (3.5b) 

Above results show that in the lower magnetic impurity band the carrier spin is almost completely 
antiparallel with the fluctuating localized spin at an M site. 

In Fig. 4, we extract the low energy part of the DOS. A magnetic impurity band forms around 
the impurity level and its DOS has similar shape to the model band. The total number of states 
in the impurity band per site is equal to x irrespective of the value of (S z ). When (S z ) = S, 
all states in the impurity band are down-spin states, whereas the impurity band is composed of 
the same number (x/2) of up- and down-spin states when (S z ) = 0. The impurity band has a 
larger bandwidth in the ferromagnetic state than in the paramagnetic state. The change in the 
bandwidth is understood as a result of the change in the effective hopping integral between Mn 
sites. When the coupling between the localized spin and the carrier spin is strong, the effective 
hopping between site i and j is reduced in proportion to cos Oij/2 where 0y is the angle between the 
localized spins at site i and j. 25 -* In a fully polarized ferromagnetic state cos#jj/2 is always unity 
while (cos 9ij/ 2) = 2/3 if the localized spins are completely randomly oriented. The result shown 
in Fig. 3 certifies that Anderson-Hasegawa theory is applicable for the magnetic impurity band in 
this case. Change in the bandwidth makes the energy of the ferromagnetic state lower than that of 
the paramagnetic state when n is small. The energy gain increases with increase of n and reaches 
a maximum at n ~ x/2. and then gradually decreases to vanishes at n ~ x as will be shown later. 
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3.3 The case of IS = -0.4 A and E M = 

Here we discuss the case of the exchange energy IS = —0.4 A and the band offset energy 
Em = 0.0. In this case the acceptor level does not appear in the dilute limit since \Em±IS\ < 0.5 A. 
To our knowledge no magnetic impurity band has been reported in A^_ x Mn x B^ l -type DMS's. 1 ' 
Therefore the present model with the parameters of \Em ± IS\ < 0.5A may correspond to II- VI- 
DMS's. Numerical results are presented in Figs. 5 ~ 7. In Fig. 5 the spin-polarized DOS is shown 
for various values of (S z ). As can be seen in Fig. 6(b), the carrier states at the M -site exist in the 
whole range of band energy. The results for P(to) and Q(oS) shown in Fig. 7 suggest that coupling 
between the carrier spin and the localized spin is not strong except at the band edges, which is 
consistent with the weak (S z ) dependence of xD m (oj) shown in Fig. 6(b). In Fig. 5 the Fermi level 
Ep for n = 0.05 is indicated by an arrow. 

3.4 Property of ferromagnetism in Em = 

In order to study the nature of carrier-induced ferromagnetism, we first calculate the Curie 
temperature T c in a very simple way. We estimate T c using the energy of the paramagnetic state 
E(0) and that of the completely ferromagnetic state E(S) as 

k B T c = ^ [ E (0) - E(S)] (3.6) 

assuming b = in Eq. (2.24). Note that the CPA condition (2.15) results in a cubic equation when 
(S'z)av = S and a quartic equation when (S z ) av = 0. 26 ~ 28 ) 

In Fig. 8 the result for T c is presented as a function of n for various values of IS/ A. We 
immediately notice that there are two different types of behavior in T c as a function of n depending 
on the size of \IS\/A. When \IS\/A is small (|75|/A < 0.3), the ferromagnetism occurs over a 
wide range of n. The Curie temperature gradually increases with the increase in n and reaches 
a broad maximum. Then it gently decreases and vanishes at a critical value n c . The maximum 
T c stays at a low value, and n c is much larger than x{= 0.05). On the other hand, when \IS\/ A 
is large (\IS\/A > 0.7), the ferromagnetism occurs in a narrow range of n (< x). The T c rises 
steeply and reaches a maximum at n x « x/2, and then it decreases rapidly. The maximum 
T c is high and n c is less than but nearly equal to x. These two different features can be seen 
clearly in Fig. 9 as well, where n c and the maximum T c are depicted as functions of n. The 
carrier density n x at which T c reaches the maximum is also shown. Two different characteristic 
features were also recognized in the phase diagrams obtained in an earlier study with Ising localized 
spins. 22 ) When \IS\/A < 0.3, the maximum T c is approximately proportional to (IS/ A) 2 . This 
suggests that in the range of \IS\/A < 0.3 the perturbative treatment on IS/ A is available and 
the RKKY-like mechanism is expected to operate for a moderate carrier density. With further 
increase in \IS\/A the maximum T c rises rapidly (0.3 < \IS\/A < 0.7), and then tends to saturate. 
For \IS\/A > 0.7, ferromagnetism is induced only when n < x and the maximum T c realizes 
at n x = x/2. The case with IS = —A studied in §3.2 belongs to this region and therefore the 
mechanism for ferromagnetism is understood from the argument above. As was shown in Fig. 4 
the width of the magnetic impurity band increases with increase of magnetization according to 
Anderson-Hasegawa mechanism. It leads to the energy gain in ferromagnetic state if the Fermi 
level lies in the impurity band. This is essentially the same mechanism with that of the double- 
exchange interaction which causes ferromagnetism in mixed valency transition metal oxides such 
as (LaCa)Mn03. 24 ' 25 ) Therefore we may say that the 'double-exchange (DE)-like' mechanism for 
ferromagnetism is operative in the impurity band when \IS\/A > 0.7. Here we point out that 
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the magnetic impurity levels appear when \IS\/A > 0.5 if Em = 0. Therefore we conclude that 
the DE-like mechanism in a magnetic impurity band becomes dominant when |/,S|/A > 0.7. In 
II- VI DMS's, the magnetic impurity level does not appear as illustrated in Fig. 5, in which the 
parameters, |/5|/A = 0.4 and Em = 0, were employed. Therefore, we may conclude that the 
DE-like mechanism is not relevant to the ferromagnetism in II- VI DMS's. 

To confirm above picture, we try to estimate T c in the strong exchange interaction limit through 
a simple argument. In this limit the impurity band has the same shape of the model density of 
states p(to) defined by Eq. (2.17) around the impurity level E a . In the completely ferromagnetic 
case the DOS of the impurity band has a finite value at the energies of E a — yfxA < u < E a + y^A 
and is given by 

D F (u) = ^xp(^SA . (3.7) 



In the paramagnetic case the DOS of the impurity band has a finite value at the energies of 

/x/2 A and is given by 



D P (u>) = 2 x Vx/2 x p - " , (3.8) 




where the factor 2 is due to up- and down-directions of the carrier spin. The result for T c based on 
the assumption is inserted in Fig. 8 as a 'LIMIT'. The maximum T c is estimated to be 

k B T c = 2(2 ~^ } v^A, (3.9) 

97T 

at n = x/2. The value of the maximum T c (= 0.0093A) for x = 0.05 is pointed by an arrow in Fig. 
9. 

It is worth noting that the Zener double-exchange mechanism for ferromagnetism is usually 
understood to be effective only when the exchange energy is larger than the width of the carrier band 
(or | JS] > 2A) in the case where magnetic ions sit on each site. In the present case the exchange 
energy is not larger than the width of the model band. Nevertheless the 'DE-like' mechanism 
functions because the magnetic impurity bandwidth is smaller than the exchange energy (or |LS| > 
v/2xA). 

3.5 Effect of the attractive potential 

We study here the role of the attractive potential in order to elucidate the origin of the carrier- 
induced ferromagnetism in III-V-based DMS's. Figures 10 ~ 12 show the results for the case where 
IS = — 0.4A and Em = — 0.6A. Although the strength of the exchange coupling is same as that 
in §3.3, an impurity band appears due to an attractive nonmagnetic local potential Em = — 0.6A. 
Furthermore, we find a strong similarity between the low energy part of DOS in Fig. 1 (IS = — A 
and Em = 0) and that in Fig. 10 (IS = — 0.4A and Em = — 0.6A). The similarity is due to 
the fact that the impurity level has the same energy, E a = — 1.25A, which is determined by the 
effective attractive potential Em + IS. When (S z ) = a magnetic impurity band with x/2 down- 
and up-spin states forms around the impurity level in both cases. Thus the impurity bands in 
both cases have a strong similarity although they do not completely agree with each other. When 
(S z ) = S, the down-spin bands agree completely with each other in both cases because the DOS's 
were determined only by the energy of antiparallel coupling Em + IS. The up-spin impurity bands 
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are completely suppressed in both cases and the states merge into the host band. Thus, for up-spin 
as well, we find a similarity in the low energy part of the DOS's shown in Figs. 1 and 10. From the 
similarity in the low energy part of the DOS, we may expect that ferromagnetism occurs through 
the same mechanism in both cases. In Fig. 13 the effect of the nonmagnetic potential Em on T c 
is presented for IS fixed to be — 0.4A. The impurity level appears for Em < — 0.1A in this case. 
When Em > 0.0 the T c stays low and n c is much larger than x, while for Em < — 0.2A high T c 
is realized and n c is less than x. In the latter region the DE-like mechanism becomes in effect. 
Figure 14 shows how the maximum T c as well as n c varies with Em and IS. The result reveals 
that the attractive nonmagnetic potential Em extends the region where the DE-like mechanism is 
effective and enhances maximum T c . The criterion for the DE-like mechanism to operate is roughly 
estimated to be IS + 0AE M < -0.6A. 

§4. Summary and Discussions 

In the present study we applied the dynamical CPA to a simple model in order to investigate 
the carrier states of DMS's in a systematic way. The present results reveal that a magnetic impurity 
band forms separate from the host band when the effective attractive potential is strong enough. 
In such a case the spin of a carrier in the impurity band couples very strongly with the localized 
spin at a magnetic impurity site. Then the hopping of the carriers among the magnetic sites aligns 
the localized spins parallel through the DE-like mechanism. The effective attractive potential 
is determined by the sum of the non-magnetic attractive potential and the exchange interaction. 
Therefore, the non-magnetic potential substantially enhances the effect of the exchange interaction. 
This mechanism works even in the case where the impurity band is not completely separate from 
the host band. This is the case in III-V-based DMS's, where the exchange interaction alone is 
not strong enough to make a pronounced impurity band. On the other hand, the absence of non- 
magnetic attractive potential in II-VI-based DMS's suppresses the magnetic impurity band and 
therefore the appearance of strong ferromagnetism. 

Here we compare the Zener double-exchange (DE) mechanism and the present DE-like mecha- 
nism that occurs in the magnetic impurity band. Zener originally proposed the DE interaction for 
(La,Ca)Mn03 where 3c?-holes hop among the magnetic ions sitting on the regular lattice sites. 34 ) 
Therefore it might be usually understood that DE mechanism is only relevant to hopping of the 
3d holes in (Ga,Mn)As. The DE mechanism for ferromagnetism is, in fact, works quite generally 
in cases where the kinetic energy (or the bandwidth) of a carrier is not greater than the exchange 
energy between a carrier spin and a localized spin. It originates from the effective reduction of the 
hopping amplitude between sites with localized spins which are not parallel with each other. The 
reduction is in proportion to cos 9/2 where 9 is the relative angle between the localized spins. 25 ' 
Only condition required for this reduction is a very strong exchange coupling between carrier spins 
and localized spins. If it is satisfied, carriers may have any character and the localized spins can 
be arranged randomly. Even localized spins are not necessary for this mechanism if strong Hund 
couplings are working between carrier spins on the same site. 35,36 ) 

In the case of III-V based DMS's, carriers are considered to have 4p character 6, 9 ) and therefore 
the bandwidth (~4 eV) is not larger than the exchange energy (IS = 0.6 ~ 0.8 eV). However the 
attractive interaction between the carrier and a Mn 2+ ion helps for the ion to bind the carrier, 
and in consequence the carrier moves in a narrow impurity band. The effective bandwidth for a 
carrier is that of the impurity band and it is smaller than the exchange energy. The DE mechanism 
becomes in effect even if the impurity band is not completely separated from the host band. We call 
this mechanism "DE-like" mechanism in order to avoid confusion with the argument assuming d 
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holes. 17 ) We note here that no Mn 3+ (d 4 configuration) states have been experimentally detected in 
(Ga,Mn) As. 3 ' 5 ' 6 ' 13, 37 ) All these experimental observations rather suggests the picture that the fixed 
valence state Mn 2+ (5 = 5/2) is realized in (Ga,Mn)As. The absence of the attractive potential 
suppresses the appearance of a magnetic impurity band and therefore the DE-like mechanism does 
not work in II-VI-based DMS's. 

The results for the n dependence of the energy gain and T c , shown in Fig. ?? are well explained 
on the basis of the DE-like mechanism. Recently Hayashi et al. reported that heat treatments 
(annealing) could vary the hole concentration while holding x constant. 15 ) They reported in 
a sample with a nominal Mn content of 5 % that T c increased from ~80K to ~95K and then 
decreased again to ~80K when rih was increased from 6.1 to 6.8 and then 7.2 xlO 19 cm" 3 . The 
carrier concentration in these cases are 2.7, 3.1 and 3.3%, respectively. The result implies that T c 
takes its maximum at nh ~ 0.6x in the sample, while the present theory predicts a maximum at 
~ 0.3x as is shown in Fig. 20(b). Study of a more realistic model taking into other effects and/or 
more detailed experimental investigations are necessary to clarify this discrepancy. 

Although we assumed the p-holes with a wide bandwidth for carriers, the present result is 
in contrast with the RKKY argument which assumes the free-carrier picture. 16 ) Infrared optical 
absorption measurements have recently revealed that the hole wave functions are nearly localized 
and the transport should be more or less hopping-like even in the most metallic Gai-^Mn^As. 38 ' 39 ) 
The nearly bound-hole picture is consistent with the present result that the p hole in the band tail 
prefers to stay on Mn sites [see Fig. ??(c)]. 

Although the exchange interaction between p hole and d spin is experimentally established to 
be anii/erromagnetic, 6 ' 12 ' 13 ) the exchange interaction was reported to be /erromagnetic in an early 
study on the polarized magnetoreflection. 40 ) In our result shown in Fig. 4, the bottom of the down- 
spin magnetic impurity band shifts to the lower energy side with increase in (S z ), whereas that of 
the host band shifts to the higher energy side. It is also worth noting that Q(u>) « —1.0 at the 
bottom of the host band (see Fig. 3 (b)). These indicate that the ferromagnetic spin coupling occurs 
at the bottom of the host band, contrary to the antiferromagnetic spin coupling in the magnetic 
impurity band. If we interpret that the optically observed band edge is not that of the impurity 
band but that of the host band, the present result may explain the experimental observation. A 
simple calculation based on the present approach consistently explains the apparent ferromagnetic 
exchange interaction experimentally observed. Details will be published elsewhere. 41 ) 

Next, we briefly mention the quantum effect of the localized spin that is ignored in the present 
work. When the calculations are done for the quantum spins with finite 5, the total number of 
states in a impurity band is not x but 2x(S + l)/(25 + 1) for parallel coupling and 2x5/(25 + 1) 
for antiparallel coupling, respectively. The result is valid as long as the band is almost empty. The 
many body effect, however, may change the situation. 42 ) Therefore, although the present result 
suggests that T r mig ht be multiplied by a factor ^2(5 + l)/(25 + 1) or ^25/(25 + 1) in the strong 
coupling limit, the effect of quantum spin is still unclear at the present stage. 

The Curie temperature obtained by using the model with randomly distributed Ising spins 22 ) 
is about three times larger than the present result. The reason may be due to the fact that the 
spin-flip effect is not taken into account in their theory. 

Throughout this work, the carriers were assumed to be degenerate. Recently the dynamical 
mean-field theory 43 ' 44 ) was applied to the present model. The result for T c is almost the same 
as the present result, which suggests that the assumption of the degenerate carriers is reasonable. 
The result will be published elsewhere together with the result of the temperature dependence of 
magnetization. 45 ) 
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The present model has succeeded in explaining the fundamental mechanism of the carrier- 
induced ferromagnetism in III-V based DMS's. There are, however, many features which exist in 
the real DMS's but are not taken into account. Those are multi-band effects, the band anisotropy, 
the spin-orbit coupling, p-d hybridization, the long-range character of the Coulomb potential, the 
direct antiferromagnetic exchange interaction between Mn ions and so on. These issues remain for 
future study. 
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Fig. 1. The result for the DMS with IS = —A and E M = 0. The density of states (DOS) as a function of lo/A for 
various values of (S z )/S. Solid line is for down-spin carrier, and dotted line is for up-spin carrier. The impurity 
levels E a = ±1.25A are indicated by dots on the abscissa. 
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The results for the DMS with 7S : 
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(1 - x)Df(u)A. (b) Af-site component of the DOS, xDf (u)A and -xD™{u)A. The thick, thin and dotted 
lines represent the cases of (S z )/S = 1.0, 0.5 and 0.0, respectively. Note the difference in the scale of vertical 
axes between (a) and (b). 
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Fig. 3. The result for the DMS with IS = — A and Em = 0. (a) Optical carrier-spin polarization P(u>) for {S z )/S = 
1.0,0.8,0.6,0.4,0.2 and 0.0. (b) Spin-coupling strength Q(u>). The thick, thin and dotted lines represent the cases 
of (S z }/S — 1.0, 0.5 and 0.0, respectively. The impurity levels E a — ±1.25A are indicated by dots on the abscissa. 



down-spin DOS 




Fig. 4. The DOS of the impurity band in the case of IS/ A = —1.0 and Em = 0. The thick, thin and dotted lines 
represent the cases of {S z )/S = 1.0, 0.5 and 0.0, respectively. The dot A indicates the impurity level for x — > 0. 
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Fig. 5. The result for IS = -0.4A and E M = 0. The DOS as a function of cu/A for various values of (S z )/S. Solid 
line is for down-spin carrier, and dotted line is for up-spin carrier. The arrow points to the Fermi level sf/A for 
n = z(=0.05). 
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Fig. 6. The result for the DMS with IS = -0.4A and E M = 0. (a) A-site component of the DOS, (1 - x)Df(w)A 
and -(1 - x)Df(cu)A. (b) M-site component of the DOS, xDf (w)A and -iD" (w)A. The thick, thin and 
dotted lines represent the cases of (S z )/S = 1.0, 0.5 and 0.0, respectively. 
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Fig. 7. The result for the DMS with IS = — 0.4A and Em = 0. (a) optical carrier-spin polarization P(u), (b) 
spin-coupling strength Q(u). 




Fig. 8. The result for Curie temperature T c /A as a function of carrier density n for various values of IS/ A with 
x — 0.05 and Em = 0. The result based on the assumption that an impurity band is similar shape of the model 
band is drawn as 'LIMIT' (see text). 
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Fig. 9. Phase diagram for Em = and x = 0.05. The critical value n c (solid line; left scale) and the maximum T c 
(dotted line; right scale) are presented as a function of \IS\/A. The carrier density n x at which T c reaches the 
maximum is included (solid line; left scale). 




Fig. 10. The result for the DMS with IS = -0.4A and E M = -0.6A. The DOS is plotted as a function of w/A for 
various values of (S z )/S. Solid line is for down-spin carrier, and dotted line is for up-spin carrier. The impurity 
level E a — —1.25 A is indicated by the dot on the abscissa. 
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Fig. 11. The result for the DMS with IS = -0.4A and E M = -0.6A: (a) ,4-site component of the DOS, (1 - 
x)Df(u)A and -(1 - x)Df(u)A, (b) M-site component of the DOS, xDf (lu)A and -xDf (u>)A. The thick, 
thin and dotted lines represent the cases of {S z )/S = 1.0, 0.5 and 0.0, respectively. 




Fig. 12. The result for the DMS with IS = — 0.4A and Em = — 0.6A: (a) optical carrier-spin polarization P(u)), (b) 
spin-coupling strength Q(w). The impurity level E a = — 1.25A is indicated by a dot on the abscissa. 



20 



Masao Takahashi and Kenn Kubo 




Fig. 13. The result for T c /A as a function of n for various values of Em /A with x = 0.05 and IS = — 0.4A. The 
result of 'LIMIT' is included (see text). 




Fig. 14. Phase diagram for various values of Em /A. The critical value n c (solid line; left scale) and the maximum 
T c (dotted line; right scale) are presented as a function of \IS\/A for various values of Em I 'A. 



